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Motivated by the anisotropic long-range nature of the interactions between cold dipolar atoms 
or molecules in an optical lattice, we study the anisotropic quantum Heisenberg model with Curie- 
Weiss- type long-range interactions. Absence of a heat bath in optical lattice experiments suggests a 
study of this model within the microcanonical ensemble. The microcanonical entropy is calculated 
analytically, and nonequivalence of microcanonical and canonical ensembles is found for a range of 
anisotropy parameters. From the shape of the entropy it follows that the Curie- Weiss Heisenberg 
model is indistinguishable from the Curie- Weiss Ising model in canonical thermodynamics, although 
their microcanonical thermodynamics differs. Qualitatively, the observed features of nonequivalent 
ensembles are expected to be relevant for long-range quantum spin systems realized in optical lattice 
experiments. 
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Cold dipolar gases have been in the very focus of ex- 
perimental and theoretical research recently [H, Q. In 
particular, dipolar gases in optical traps have been sug- 
gested as laboratory realizations of lattice spin models 
where the coupling parameters can be tuned freely, allow- 
ing for the realization of many Hamiltonians of interest 
in condensed matter physics [3[. 

After switching off the cooling in such an experiment, 
total energy and number of atoms are conserved to a very 
good degree. As a consequence, a statistical description 
of such a lattice spin model should make use of the micro- 
canonical ensemble. For systems with short-range inter- 
actions, the choice of the statistical ensemble is typically 
of minor importance and could be considered a finite- 
size effect: differences between, say, microcanonical and 
canonical expectation values are known to vanish in the 
thermodynamic limit of large system size, and the var- 
ious statistical ensembles become equivalent [4[. In the 
presence of long-range interactions this is in general not 
the case, and microcanonical and canonical approaches 
can lead to different thermodynamic properties even in 
the infinite-system limit [5j. In the astrophysical con- 
text, nonequivalence of ensembles and the importance 
of microcanonical calculations have long been known for 
gravitational systems 0, 0] • 

In condensed matter physics, most systems are coupled 
to an environment, and therefore the canonical or grand- 
canonical ensembles are the ones that appropriately de- 
scribe the experimental situation of interest. Moreover, 
screening effects lead in general to interactions which are 
effectively of short range, and hence equivalence of en- 
sembles usually can be taken for granted. As a conse- 
quence, calculations of thermodynamic quantities can be 
done in the ensemble that is the most convenient one, 
which appears to be the canonical or grandcanonical, but 
never the microcanonical, one. Owing to these facts, lit- 



tle is known about such systems in the microcanonical 
ensemble. Only recently, a number of toy models, con- 
sisting of long-range coupled classical spin variables, has 
been studied (see [8| for a review). The study of these 
strongly simplified but analytically solvable models has 
been very fruitful towards the aim of understanding gen- 
eral dynamical and thermodynamical properties of clas- 
sical systems with long-range interactions. 

Much less is known about the peculiarities of quantum 
spin systems with long-range interactions, and, in partic- 
ular, about equivalence or nonequivalence of ensembles in 
this context. It is the aim of this Letter to contribute to- 
wards the understanding of such systems, with a focus 
onto the microcanonical setting as encountered in exper- 
iments with dipolar gases in optical lattices. 

To this purpose, we study the anisotropic quantum 
Heisenberg model with Curie- Weiss-type long-range in- 
teractions in the microcanonical ensemble. Such Curie- 
Weiss- type interactions, where each spin is interacting 
with every other at equal strength, are clearly an ideal- 
ization of the actual interactions of dipolar atoms which 
decay like r~ 3 with the interparticle distance r. However, 
it is known that Curie- Weiss-type models faithfully re- 
produce many properties of algebraically decaying long- 
range interactions qualitatively, and to some extent even 
quantitatively foL flo|. 

In this Letter the result of an exact, analytic calcula- 
tion of the thermodynamic limit of the microcanonical 
entropy of the anisotropic quantum Heisenberg model 
with Curie- Weiss-type interactions is reported. Depend- 
ing on the choice of the anisotropy parameters in the 
Hamiltonian operator, a concave entropy function is 
found in some cases, and a nonconcave one in others. 
Correspondingly, equivalence of the microcanonical and 
the canonical ensemble holds in the first case, but not in 
the second. 
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The relevance of the reported results is twofold. First, 
the observation of nonequivalent ensembles in long-range 
quantum spin systems demonstrates that, under the ex- 
perimental conditions realized in cold dipolar gases in 
optical traps, the choice of the statistical ensemble is of 
paramount importance. Consequently, a statistical in- 
terpretation of the results of such experiments has to 
go beyond usual canonical thermodynamics. In partic- 
ular, differences between microcanonical and canonical 
expectation values do not diminish in importance with 
increasing system size. This is in sharp contrast to mi- 
crocanonical computations for ideal Bose gases in traps 
[ll|], where equivalence of ensembles holds in the ther- 
modynamic limit. Second, the reported calculation also 
illustrates that cold dipolar gases in optical traps are ex- 
cellent laboratory systems in which long-range effects like 
the nonequivalence of statistical ensembles or the nega- 
tivity of microcanonical response functions can possibly 
be tested. Since these effects occur only for a certain 
range of values of the anisotropy parameters, it is of par- 
ticular importance that coupling constants (and therefore 
anisotropy parameters) in cold atom experiments can be 
tuned with a high level of control, rendering such systems 
an ideal laboratory for the study of these fundamental is- 
sues of thermostatistical physics. 

Anisotropic quantum Heisenberg model — The model 
consists of N spin- 1/2 degrees of freedom, each of which 
is interacting with every other at equal strength. The 
corresponding Hilbert space % = ((D 2 ) 0Ar is the tensor 
product of N copies of the spin- 1/2 Hilbert space (D 2 , 
and the Hamilton operator is given by 

I N N 

k,l=l k=l 

(i) 

The a% are operators on H and act like the a component 
of the Pauli spin- 1/2 operator on the fcth factor of the 
tensor product space H, and like identity operators on all 
the other factors. The resulting commutation relation is 

[at, of] = 2i 5 k ,i e aM al, a, /3 G {1, 2, 3}, (2) 

where 5 denotes Kronecker's symbol and e is the Levi- 
Civita symbol, h is the strength of an external magnetic 
field orientated along the 3 axis, and the constants Ai, A2, 
and A3 determine the coupling strengths in the various 
spatial directions and allow us to adjust the degree of 
anisotropy. Note that it is explicitly shown in [3] that 
anisotropic quantum Heisenberg models are among the 
systems that can be engineered with cold polar molecules 
in optical lattices. 

Special choices for the coupling constants in (pQ) yield, 
for example, (a) the isotropic Heisenberg model, Ai = 
A2 = A3, (b) the Ising model, Ai = = A2, (c) the 
isotropic Lipkin-Meshkov-Glick model, Ai = A2 and A3 = 
0. For these special cases, the Hamiltonian (pQ) can be ex- 
pressed in terms of S 2 and 63, i.e., the square and the 3 



component of a collective spin operator S = ^cr fc /2. As 
a consequence, an angular momentum eigenbasis simul- 
taneously diagonalizes and S3 and the model can be 
solved by elementary means. 

Here we consider the coupling constants Ai, A2, and 
A3 to be nonnegative, but otherwise arbitrary, and in 
this case the model is known to display a transition from 
a ferromagnetic to a paramagnetic phase in the canonical 
ensemble. The exact expression for the canonical Gibbs 
free energy g as a function of the inverse temperature 
/3 = 1/T [12] and the magnetic field h is known for this 
model (and, in fact, for a larger class of systems) and can 
be found, for example, in [13|. 

Microcanonical entropy. — In thermodynamics, the en- 
ergy e is the variable conjugate to the inverse tempera- 
ture /?, and the magnetization m is conjugate to —j3h. So 
in the same way that g(/3, h) represents the fundamental 
quantity of the quantum Heisenberg model in the canon- 
ical ensemble, the microcanonical entropy s(e, m) serves 
as a starting point for a microcanonical description in 
the thermodynamic limit. However, for a pair of vari- 
ables (e, m) corresponding to the pair of noncommuting 
operators (Hq,M = 2 S3), it is not even well established 
how to define a quantum microcanonical entropy, sym- 
bolically given by 

s N (e,m) = ^]nTr[6(Ne-H )5(Nm-M)]. (3) 

Note that the symbolic expressions make little mathe- 
matical sense and require some physically reasonable reg- 
ularization. Extending a suggestion by Truong [14[ to 
interacting systems, the definition 

s N (e,m) = ^\n^Tr[P Ho (e)P M (fh)} 5 A (e-e)5 A (fh-m) 

e,fh 

(4) 

seems to be physically reasonable, but difficult to apply 
in practice. Here, e and fa denote eigenvalues of the op- 
erators Ho/N and M/N, respectively. 6 a is the charac- 
teristic function of the interval [— A,0], i.e., 5&(x) — 1 if 
x G [— A,0], and zero otherwise. P# (e), Pm(™>) denote 
the eigenprojections of the operators Ho and M belong- 
ing to the eigenvalues e and m, respectively. 

We here report results for s(e, m) = liniAr^oo Sjv(e, m) 
obtained by using a different regularization. The ana- 
lytic calculation of s uses, among others, some ingredi- 
ents from a related canonical calculation by Tindemans 
and Capel [15j. Details will be reported elsewhere, but 
the main steps of the calculation can be sketched as fol- 
lows, (i) The deltas in ((3]) are replaced by their Fourier 
integral representations, (ii) The Lie- Trotter formula 
is applied to separate the resulting exponential of the 
Hamiltonian into exponentials of the type exp{c Q! (5' a ) 2 }, 
a G {1,2,3}, with constants c a and collective spin com- 
ponents S a . (hi) These exponentials are transformed 
into exponentials exjp{c a S a } by applying the Hubbard- 
Stratonovich trick. The trade-off for steps (ii) and (iii) is 
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a 3n+2-dimensional integral, to be considered in the limit 
n — >• oo. The advantage, however, is that the Hilbert 
space trace of exp{c a S a } factorizes into traces over the 
single-spin Hilbert spaces (D 2 , which can be easily per- 
formed, (iv) The resulting high-dimensional complex in- 
tegral can be solved in the thermodynamic limit N — >• oo, 
for example, by the method of steepest descent. 

The final result for the microcanonical entropy of the 
Curie- Weiss anisotropic quantum Heisenberg model in 
the thermodynamic limit is 

s(e,m) =ln2 - ^[1 - /(e,ra)]ln[l - /(e,m)] 

1 2 (5) 
--[l + /(e,m)]In[l + /(e,m)] 

with 

,(..„,) = ^(l -£)-£. (6) 

and A^ = max{Ai,A2} [l6j], where s(e, m) is defined on 
the subset of R 2 for which 

< m 2 (X_ L - A 3 ) -2e < X± and 2e < -m 2 A 3 . (7) 

The result is remarkably simple, in the sense that an 
explicit expression for s(e, m) can be given. This is in 
contrast to the canonical ensemble, where g(/3, h) is given 
implicitly as the solution of a maximization [13|. Plots 
of the domains and graphs of s(e, m) are shown in Fig. [T] 
for a number of coupling strengths Aj_, A3. 

Nonequivalence of ensembles. — On a thermodynamic 
level, equivalence or nonequivalence of the microcanon- 
ical and the canonical ensembles is related to the con- 
cavity or nonconcavity of the microcanonical entropy [H| . 
By inspection of rows three to seven in Fig. [T] [or by sim- 
ple analysis of the results in ©-©], the entropy s for 
A_l > A3 is seen to be a concave function on a domain 
which is a convex set. For A^ < A3, the domain is not a 
convex set and therefore the entropy is neither convex nor 
concave. In the latter case, microcanonical and canoni- 
cal ensembles are not equivalent, in the sense that it is 
impossible to obtain the microcanonical entropy s(e,m) 
from the canonical Gibbs free energy g(/3,h), although 
the converse is always possible by means of a Legendre- 
Fenchel transform. 

The physical interpretation of ensemble equivalence is 
that every thermodynamic equilibrium state of the sys- 
tem that can be probed by fixing certain values for e 
and m can also be probed by fixing the corresponding 
values of the inverse temperature /3(e, m) and the mag- 
netic field h(e,m). In the situation A^ < A3 where 
nonequivalence holds, this is not the case: only equilib- 
rium states corresponding to values of (e, m) for which 
s coincides with its concave envelope can be probed by 
fixing macrostates corresponding to other values 




FIG. 1. (Color online) Domains (left) and graphs (right) of 
the microcanonical entropy s(e,m) of the anisotropic quan- 
tum Heisenberg model for some combinations of the couplings 
A_l,A 3 . From top to bottom: (A_l,A 3 ) = (1/4,1), (9/10,1), 
(1,1), (1,9/10), (1,1/2), (1,1/5), (1,0). For the domains, the 
abscissa is the energy e and the ordinate is the magnetization 
m, and the entropy is defined on the shaded area. 
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of (e, m), however, are not accessible as thermodynamic 
equilibrium states when controlling temperature and field 
in the canonical ensemble. In this sense, microcanonical 
thermodynamics can be considered not only as differ- 
ent from its canonical counterpart, but also as richer, 
allowing to probe equilibrium states of matter which are 
otherwise inaccessible. The realization of a long-range 
quantum spin system by means of a cold dipolar gas in 
an optical lattice offers the unique and exciting possibil- 
ity to study such states in a fully controlled laboratory 
setting [l7|. 

Thermodynamic equivalence of models. — Let us leave 
aside for a moment the question of experimental realiza- 
tion and discuss a different kind of equivalence specific to 
the anisotropic quantum Heisenberg model. It had been 
observed already in the 1970s that the isotropic Heisen- 
berg model and the Ising model are thermodynamically 
equivalent in the sense that their canonical free energies 
coincide [18| . One can verify by Legendre-Fenchel trans- 
forming the entropy in (j5]) that the same is in fact true for 
all coupling strengths satisfying ^ A3. Geometrically, 
this thermodynamic equivalence corresponds to the fact 
that the entropies s(e, m) for those couplings share the 
same concave hull (which is equal to the entropy of the 
isotropic Heisenberg model plotted in row three of Fig. H}. 
Identical concave hulls of entropies imply, however, iden- 
tical canonical free energies, and hence thermodynamic 
equivalence of the two models follows in the canonical 
ensemble. Remarkably, however, thermodynamic equiv- 
alence does not hold in the microcanonical ensemble, as 
is obvious from the different shapes of entropies in rows 
one to three of Fig. HJ 

Discussion. — The microcanonical entropies (j3j) and 
(|4]) discussed in this Letter describe the physical situ- 
ation of fixed energy e and magnetization m. In a cold 
atom experiment, energy is conserved to a very good de- 
gree due to the absence of a heat bath. For apolar gases 
where s-wave scattering is dominant, the total magneti- 
zation is also fixed, and the resulting short-range inter- 
acting microcanonical spin systems have been discussed 
in [19[. For dipolar gases where long-range interactions 
are present and nonequivalent ensembles can occur, the 
magnetization is not conserved in general (unless an ex- 
perimentalist comes up with an ingenious trick the author 
is not aware of). More easily, nonequivalence of ensem- 
bles could be observed in long-range quantum spin sys- 
tems undergoing a temperature-driven first-order transi- 
tion. In this case, nonequivalence is signalled by a non- 
concave microcanonical entropy s(e), corresponding to 
conservation of energy e, but fluctuating magnetization. 

Although the anisotropic quantum Heisenberg model 
discussed in this Letter is among the systems which can 
be engineered with cold polar molecules in optical lattices 
Q, this model is chosen here not for its particular fea- 
tures, but to illustrate general, and possibly even generic, 
properties of long-range interacting quantum spin sys- 



tems: nonconcave entropies, nonequivalence of statisti- 
cal ensembles, or other phenomena like negative micro- 
canonical response functions must be expected to show 
up under the experimental conditions realized in exper- 
iments with cold dipolar gases in optical lattices in gen- 
eral. Finally, note that dipolar atoms or molecules are 
not the only possible realization of long-range quantum 
spin systems in optical lattices: Following a suggestion by 
O'Dell et al. [20], long-range interactions decaying with 
the interparticle distance r like r _1 can be engineered by 
shining appropriately tuned laser light onto atoms, even 
in the absence of a permanent dipole moment. 

Summary. — A calculation of the microcanonical en- 
tropy of the anisotropic Curie- Weiss quantum Heisenberg 
model was reported. The results illustrate peculiarities 
of long-range quantum spin systems, like nonconcave en- 
tropies and nonequivalence of statistical ensembles. The 
microcanonical setting models the conditions relevant for 
experiments with dipolar gases in optical lattices. The 
results point out the importance of nonstandard thermo- 
dynamics beyond the canonical ensemble for such exper- 
iments on the one hand, and on the other hand suggest 
the use of optical lattice experiments for the study of 
fundamental issues of thermostatistics. 
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